Spin polarization and g-factor enhancement in graphene nanoribbons in magnetic field 
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We provide a systematic quantitative description of spin polarization in armchair and zigzag 
graphene nanoribbons in a perpendicular magnetic field. We first address spinless electrons within 
the Hartree approximation studying the evolution of the magnetoband structure and formation of 
the compressible strips. We discuss the potential profile and the density distribution near the edges 
and the difference and similarities between armchair and zigzag edges. Accounting for the Zeeman 
interaction and describing the spin effects via the Hubbard term we study the spin-resolved subband 
structure and relate the spin polarization of the system at hand to the formation of the compressible 
strips for the case of spinless electrons. At high magnetic field the calculated effective g-factor varies 
around a value of (<?*} ~ 2.25 for armchair nanoribbons and (<?*) ~ 3 for zigzag nanoribbons. An 
important finding is that in zigzag nanoribbons the zero-energy mode remains pinned to the Fermi- 
energy and becomes fully spin-polarized for all magnetic fields, which, in turn, leads to a strong spin 
polarization of the electron density near the zigzag edge. 
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I. INTRODUCTION 

Investigation of effects of electron interaction and spin 
in graphene at high magnetic field represents one of the 
frontiers in the graphene research. Even though many 
aspects of the magnetoconductance of graphene related 
to the formation of unconventional Landau level spectra 
and the anomalous Hall effect are well understood theo- 
retically and confirmed experimentally^—, there are still 
a number of questions awaiting their resolution. One of 
these questions which is extensively debated in the cur- 
rent literature is the origin of the splitting of the lowest 
Landau level and the emerging of insulating state at the 
Dirac point Even though the precise origin of this 
state is under current debate, it is generally believed 
that it is related to electron-electron interaction and 
spin effects. The importance of electron interaction was 
also outlined for higher Landau levels^. Recently, spin- 
splitting in graphene and bilayer graphene in high mag- 
netic field was experimentally analyzed by Kurganona et 
alJ£, who found that g-factor in graphene is enhanced, 
and attributed this to electron-electron interaction ef- 
fects. The spin-splitting of the states in graphene 1 ^ and 
graphene quantum dots 1 ^ was also studied in a parallel 
magnetic filed. 

Motivated by this interest to the electron interacton 
and spin effects in graphene in the high magnetic field, 
in the present paper we study the spin polarization and 
enhancement of the g-factor in graphene nanoribbons 
(GNRs). Note that various aspects of electron and spin 
interactions in high magnetic have been extensively stud- 
ied in conventional semiconducting quantum wires de- 
fined in two dimensional electron gas (2DEG)i2~— . One 
of the motivations for such studies is related to advances 
in semiconductor spintronics utilizing the spin degree 



of freedom for adding new functionalities to electronic 
devices^ Some of proposed and investigated devices for 
spintronics and quantum computation applications op- 
erates in the edge state regime ; 23 : 24 which obviously re- 
quires a detailed knowledge of the structure of the states 
in a quantum wire or at the edge of the 2DEG. The prop- 
erties and detailed information about propagating states 
at the boundaries are also essential for interpretation 
of experiments in various electron interferometers in the 
quantum Hall regime 2 ^ - — . Because graphene represents 
a very promising system for implementation of many de- 
vices and concepts for spintronics and quantum infor- 
mation processing applications utilizing the edge state 
transport regime a detailed knowledge of the density and 
potential profiles near the edges as well as spin prop- 
erties are important for understanding and designing of 
such devices. 



The paper is organized as follows. In Sec. II we present 
a formulation of the problem, define the Hamiltonian and 
briefly outline the self-consistent computational scheme. 
The results and discussion are presented in Sec. II. Sec- 
tion IIA discusses the potential profile and the charge ac- 
cumulation near the edges in ribbons of different widths 
and edge terminations. Section IIB is devoted to the case 
of spinless electrons focussing on the formation of com- 
pressible strips and evolution of the magnetoband struc- 
ture. Finally, based on the results of Sec. IIB, Sec. IIC 
discusses the spin splitting and the enhancement of the 
g-factor for the case of electrons with spin. The conclu- 
sions of the work are presented in Sec. III. 
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II. MODEL 

We consider an infinite GNR of the width W, located 
in an insulating substrate with the relative permittivity 
e r and subjected to the perpendicular magnetic field B, 
see inset to Fig. (Ub). A metallic back gate situated 
at the distance d from the ribbon is used to tune the 
Fermi energy in order to change an electron concentra- 
tion in the GNR. The system is modeled by the p-orbital 
tight-binding Hubbard-type Hamiltonian in the mean- 
field approximation, H = H' + H^, which is shown to 
describe carbon electron systems in good agreement with 
the first-principles calculation! ] 28 ' 29 , 



H a = - ^ tr.r+Aata a r+A,, 
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where a, a' correspond to two opposite spin states 'f, 
J,; the summation runs over all sites r = (x, y) of the 
graphene lattice, A includes the nearest neighbors only. 
The magnetic field is included in a standard way via 
the Pierel's substitution, t r ,r+A = to cxp(i27r0 r!r+ A/</ > o) 
where (f> r . r +A = J^ +A A • dl, with A being the vector 
potential, and (f>o = h/e being the magnetic flux quan- 
tum, to = 2.7 eV. (In our calculations we use the Landau 
gauge, A = (— By, 0)). The first two terms in Eq. ([T]) 
correspond to the non-interacting part of the Hamilto- 
nian, with the first term describing the kinetic energy 
of the electrons on a graphene lattice. The second term 
describes the Zeeman energy triggering the spin-splitting 
in the magnetic field, V% = ±^gpbB, where +(— ) signs 
correspond to the opposite spin states t (4-); 9 = 2 is the 
bare g-factor of pristine graphene, and the Bohr magne- 
ton fib = eh/2m e . The two last terms in Eq. ([T]) describe 
the electron interaction. The long-range Coulomb inter- 
action between induced charges in the GNR is given by 
the standard Hartree term, 
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where n r = n r ^ + n r j, is the total electron density, and 
the second term corresponds to a contribution from the 
mirror charges. The last term in the Hamiltonian ([1]) 
corresponds to the Hubbard energy, 



V5 (r) = U[n 



27' 



(3) 



and describes repulsion between electrons of the opposite 
spins on the same site. The number of excess electron at 
the site r reads. 



/oo 
p a (r,E)f FD (E,E F )dE-m, 
-oo 



(4) 



where p a (r,E) = -iS[G ff (r,r,£)] is the energy depen- 
dent local density of states (LDOS) at rero temperature; 
G a (r, r, E) is the Green's function in the real space repre- 
sentation of an electron of the spin a residing on the site 
r, Ep = eV g is the Fermi energy which value is adjusted 
by the gate voltage, and n ions = 1/Aq = 3.8 x 10 19 m -2 
is the positive charge background of ions (Aq = ^pa 2 , is 
the area per one C atom and ao = 0.142 nm is the C-C 
distance). Equations J!])-© are solved self-consistently 
using the Green's function technique in order to calculate 
the band structure, the charge density and the potential 
distribution)^ - — For a given potential distribution we 
compute the conductance using the Landauer formula 



G°(E F ) = e - J T°{E) 



df FD {E-E F ) 
dE 



dE, (5) 



where T a (E) is the total transmission coefficient for elec- 
trons with spin cr, and f F n is the Fermi-Dirac distribu- 
tion. 



III. RESULTS AND DISCUSSION 

A. Potential profile and the charge accumulation 
near the edges in ribbons of different widths and 
edge terminations. 

In the present study we aim at the description of 
spin polarization in realistically wide GNRs. In current 
experiments^— the widths of nanoribbons are W ~ 100 
nm and the magnetic field reaches B < 60 T, correspond- 
ing to the ratio W/Ib ~ 10, with Ib = y/h/eB being the 
magnetic length. At the same time, because of computa- 
tional limitations, it is difficult to treat ribbons of widths 
exceeding 50 nm. Therefore in our calculations we re- 
scale the system using ribbons of smaller width W ~ 30 
nm subjected to higher fields (up to B ~ 350 T) in order 
to keep the ratio W/Ib ~ 10 in accordance with typi- 
cal experiments)^— In very narrow ribbons the quan- 
tum confinement effects can dominate ribbon's electronic 
properties. Thus, a concern might arise whether the ob- 
tained results remain valid for realistically wide ribbons. 
In the present section we investigate how nanoribbon's 
electronic properties such as the density distribution and 
the potential profile in the vicinity of the edge change 
with the increase of the ribbon's width, and find that the 
width of W ~ 30 nm is already sufficient to capture all 
essential features of a wide ribbon or even a semi-infinite 
graphene sheet. 

In our study we consider both types of edges, armchair 
and zigzag. We will demonstrate in the subsequent sec- 
tions that main features in the spin polarization of the 
electron density and the enhancement of the <?-factor are 
rather similar for both types of edges. There is how- 
ever an important difference between them which can be 
traced to the presence of the zero-energy mode (ZEM) 
residing at the edge of the zigzag GNRs^i In the present 
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FIG. 1. (Color online). (a),(c) The self-consistent charge den- 
sities and (b),(d) potentials for the armchair and zigzag rib- 
bons of widths H / =10,30 and 50 nm calculated in the Hartree 
approximation at zero magnetic field. Only half of the ribbon 
is shown. The insets in (a),(c) show the dispersion relation for 
nanoribbons of different widths W=10, 30 and 50 nm from 
the left to the right for several lowest subbands.The inset in 
(b) shows a schematically geometry of the device; d = 30 nm. 
E F =5 t. t = 2.7 eV. a = 0.142 nm. 



is noteworthy that the electron density distribution ob- 
tained from the electrostatics (i.e. due to the Hartree 
potential Vfj(r), Eq. <j2j) ) is not altered significantly by 
magnetic field^. The charge densities and potentials stay 
qualitatively the same as the nanoribbon width increases 
and exhibit practically no difference for 30- and 50-nm 
wide nanoribbons. We thus conclude that the transverse 
confinement docs not change substantially for nanorib- 
bons wider than ~ 30 nm, and the width W = 30 nm 
is sufficient to describe realistically wide ribbons or even 
an edge of a graphene sheet. 

Let us now focus on a difference in the potential pro- 
files and the electron density distributions in a vicin- 
ity of a ribbon edge for armchair and zigzag ribbons. 
Both ribbons show strong electron accumulation near 
the edges but this accumulation is stronger in the zigzag 
GNRs. The corresponding potential profiles for armchair 
and zigzag ribbons have different shapes near the edges. 
For the armchair ribbon, the potential has a triangular 
shape, see Fig. [TJb). This was predicted and explored 
previouslyj 39 i 40 The triangular shape of the potential is 
related to the hard-wall confinement. It is noteworthy 
that a similar triangular shape of a potential is exhibited 
by eleaved-edge overgrown quantum wires where elec- 
trons also experience a hard- wall confinement 1 20 i 41 . 

The potential profile for the case of zigzag ribbon ex- 
hibits somehow different features. As in the case of 
the armchair GNRs, it gradually decreases towards the 
boundaries to form a well in the vicinity of the edges. 
However, in the close proximity to the edges, it raises 
up and crosses the Fermi energy, see Fig. [U^d). We re- 
late this feature to the zero-energy mode (ZEM) that 
traps charges. The zero-energy mode is manifested as 
disperseless energy level pinned to Ep in the ranges 
k e [§, l] and k E [§, l] • It is these trapped 

charges that raise the potential at the edges. They ef- 
fectively repulse excess charges induced in the ribbon by 
the gate and prevent the triangular well to form near 
the boundary. Therefore, the difference in the charge ac- 
cumulation and potential profiles in the armchair and 
zigzag ribbons occurs due to topological property of 
the zigzag edge termination supporting the zero-energy 
mode. It is important to stress that this difference per- 
sists into the high-field regime and can not be addressed 
by semi-classical approaches like in Ref. |39|. 



section we will demonstrate that the ZEM leads to a dif- 
ferent features in the potential and charge density pro- 
files near the edges for the cases of armchair and zigzag 
ribbons. 

Figure Q] shows the self-consistent charge density distri- 
butions and potential profiles for the armchair and zigzag 
nanoribbons of various widths W=10, 30, 50 nm. In all 
calculation the distance between the GNR and the gate 
is d = 30 nm. The calculations are performed in the 
Hartree approximation for spinless electrons at zero field 
(i.e. Vy and V% are set to in the Hamiltonian (JTJ). It 



B. LDOS, magnetobandstructure and formation of 
compressible strips for spinless electrons 

Spin polarization in conventional quantum wires is re- 
lated to the formation of compressible strips^ in the case 
of interacting spinless electrons *£r— In this section we 
therefore outline the electronic and transport properties 
of armchair and zigzag nanoribbons in the Hartree ap- 
proximation for spinless electrons (i.e. disregarding the 
Hubbard and the Zeeman interactions, Vz = Vu = 0, in 
the Hamiltonian Eq. (JTJ) focussing on the formation of 
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FIG. 2. (Color online). (a),(d) The LDOS and (b),(e)the DOS 
at the Fermi energy, and (c),(f) the two-terminal conductance 
calculated in the Hartree approximation for the armchair and 
zigzag ribbons (left and right panels respectively). The LDOS 
is shown for one graphene sublattice, say A; The LDOS for the 
sublattice B is symmetric with respect to the ribbon's axis. 
The regions with higher LDOS point to higher probability 
to find an electron there and correspond to positions of the 
compressible strip. N c and -/V e in (c) and (f) mark the number 
of occupied electron subbands in the center and near the edge 
of the ribbon respectively. The width of ribbons is W = 30 
nm corresponding to 242 and 141 carbons in the cross section 
of armchair and zigzag ribbons, respectively. Ef=5 t. The 
temperature T = 10 K. 



the compressible strips. 

Figures [2ja) (b) show the local density of states 
(LDOS), LDOS(y) = - J p{y,E) dfFD f E EF) dE, and 
the total density of states (DOS), DOS = J2 y LDOS(y), 
at the Fermi energy as a function of the magnetic field 
for the armchair and zigzag ribbons. (Note that LDOS is 
shown for one sublattice only.) It is noteworthy that the 
DOS can be accessible via magneto-capacitance or mag- 
netoresistance measurements similarly to conventional 
semiconductor structures defined in 2DE G 42 ' 43 . The 
structure of the LDOS and DOS can be understood from 
an analysis of the magnctosubband structure. We outline 
below the main features of the subband structure for the 
armchair and zigzag ribbons focussing on the differences 
and similarities between them as well as on formation of 
the compressible strips in the ribbons. (Note that evo- 
lution of the band structure for the case of the armchair 
GNRs was discussed by Shylau et al^). 



Left panels of Figs. [3] (a)-(c) show the band structure 
of armchair graphene nanoribbons for spinless Hartree 
electrons for three representative magnetic fields. Flat 
regions in the band diagrams correspond to the Landau 
levels in bulk graphene, and dispersiveness states close to 
the GNRs boundaries represent edge states correspond- 
ing to classical skipping orbits. Figure [3] (a) shows the 
band diagram for a magnetic field B = 140 T when the 
two lowest Landau levels, LL0 and LL1, arc filled. The 
first Landau level LL1 is pinned to the Fermi energy thus 
forming a compressible strip in the center of the GNR. 
The strip is called compressible when the electron density 
can be easily redistributed in order to effectively screen 
the external potential. We define a compressible strip as 
a region where the dispersion lies within the energy win- 
dow \E — Ep\ < 27rfcg r 18 ' 19 ' 31 i 44 because in this energy 
window the states are partially filled i.e., < Jfd < 1 
and thus the electron density can be easily changed. 

Because a graphene ribbon has abrupt edges the self- 
consistent potential forms the triangular wells near edges 
as discussed in the previous section. As a result the cen- 
ter of the ribbon and its edges depopulate in magnetic 
field differently. Namely, as the magnetic field increases, 
the subbands first depopulate in the center and then near 
the edges. For example, the second subband (i.e. LL1) 
is pinned to the Fermi level in the ribbon center in the 
interval B « 60 — 170 T. Within this interval it forms 
the compressible strip which is manifested itself as the 
enhanced LDOS and DOS in Figs. HJa) and (b). How- 
ever, the LL1 stays populated near the edges in a wider 
magnetic field range, B w 80 — 290 T . In the field in- 
terval B w 170 — 230 T the LL1 is depopulated in the 
center (see Fig. [3](b)). As a result the LDOS and DOS 
are practically zero, see Figs. d](a),(b). When the mag- 
netic field increases to B « 230 T, the lowest Landau 
level, LL0, is pushed up in energy and gets pinned to 
the Fermi energy. This again leads to a formation of the 
compressible strip in the middle of the wire (see Fig. [3J 
(b)) and to the enhancement of the LDOS and DOS at 
the Fermi energy as seen in Figs. [2] (a),(b). At this B, 
two different LLs are at Ep and contribute to electron 
transport: LL0 in the center and LL1 near the edges. 

The graphene ribbons with the zigzag edge termination 
exhibits features of the magnctosubband structure simi- 
lar to those for the armchair terminated ribbons. This is 
illustrated in Figs. [3J (d)-(f) which show the band struc- 
ture of the zigzag GNRs for three representative mag- 
netic fields, B = 90 T, 220 T and 300 T. As for the 
armchair GNRs, these fields correspond respectively to 
the cases when the LL1 is pinned to Ep in the middle of 
the GNR (Fig. [3J (d)), LL1 is depopulated in the middle 
of the GNR (Fig. [3] (e)), and LL0 is pinned to Ep in 
the middle of the GNR (Fig. [3] (f)). However, there are 
several striking differences between the zigzag and arm- 
chair ribbons manifested in their LDOS, DOS and the 
subband structure. First, strong electron accumulation 
and formation of the compressible strip takes place near 
the very ribbon's boundaries over all the range of mag- 
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FIG. 3. (Color online). The dispersion relation for (a)-(c) the armchair and (d)-(f) the zigzag ribbons. Left and right panels 
in every plot correspond respectively to the Hartree and Hubbard approximations. Because of the symmetry, only a half of 
the band diagrams is shown. The red and blue dotted lines show respectively spin-up and spin-down components. Vertical 
dot-dashed lines in (d)-(f) mark K and K' points of the first Brillouin zone of graphene in the zigzag ribbon. In armchair 
ribbon, K and K' points coincide at k — 0. In (d)-(f), "ZEM" marks the zero-energy mode that stays pinned to the Fermi 
energy. 



netic fields studied, see enhanced LDOS at y ~ in Fig. 
|2Jd). Because of this the total DOS in the zigzag GNR 
never drops to zero over all the range of magnetic fields, 
see Fig [He). (Note that Fig. [2{d) shows the LDOS for 
the sublattice A which is enhanced at the one edge of the 
ribbon. The sublattice B has the enhanced LDOS near 
the opposite edge of the ribbon). 

Inspection of the magnetoband structure reveals that 
these features are caused by the ZEM in the zigzag GNR 
discussed in the previous section. This mode (marked 
by "ZEM" in Figs. [2(d)-(f)) always stays pinned to the 
Fermi energy because of the high density of states for 
electrons that this mode accommodates. It is important 
to stress that the pinning of this mode to E = Ep is 
the result of the electron-electron interaction, and the 
pinning effect is apparently absent in the one-electron 
description where this mode the is always situated at 
E = 022. 

Figures [Hc),(f) show evolution of the two-terminal 
conductance as the magnetic field increases for the arm- 
chair and zigzag ribbons. In contrast to the conven- 



tional semiconductor quantum wires and quantum point 
contacts exhibiting a step-like conductance^, the con- 
ductance of GNRs reveals non-monotonic decrease with 
bumps coexisting with quantized plateau regions of mul- 

9 2 

tiples =£-. Note that Q changes by two conductance 
quantum between plateaus due to valley degeneracy of 
graphene. The origin of the bumps in the conductance 
was discussed by Shylau et al^- and was related to the 
interaction-induced modifications of the band structure 
leading to the formation of compressible strips in the mid- 
dle of GNRs. 



C. Spin splitting and the enhancement of the 
(/-factor 

Let us now turn to the case of electrons with spin 
and analyze how the Hubbard and the Zeeman interac- 
tions modify the magnetosubband structure of graphene 
nanoribbons leading to the spin polarization and en- 
hancement of the g-factor. For the Hubbard constant 



G 



we choose U = t which corresponds to the estimation of 
Refsi 29 ' 46 . Note that the recent work2£ predicts a some- 
how larger value, U ~ 3t. While the utilization of larger 
value of U leads to some quantitative differences with 
the results presented below, they remain qualitatively the 
same and the conclusions are not affected. We also stress 
that while the discussion below is focused on the case of 
high magnetic field when the two lowest Landau levels, 
LL1 and LLO, are occupied, the similar conclusion remain 
valid for lower fields as well. 
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FIG. 4. (a),(g) The local spin polarization of the charge den- 
sity, (b),(h) the unit cell occupancy, (c), (i) the spin polar- 
ization P, (d),(j) the DOS, (e), (k) the effective g-factor, and 
(f),(l) the conductance as a function of magnetic field calcu- 
lated in the Hubbard approximation for the armchair (a)-(f) 
and zigzag (g)-(l) ribbons. The red lines with filled dots show 
spin-up component and blue line with open squares is for spin- 
down. Arrows in (f) and (1) mark magnetic fields for the band 
structures shown in Fig. [3] 



Let us start with the case of the armchair ribbons. 
Figures H] (a)-(c) show the local spin polarization of 

the charge density, ((y) = ^fa]+n4faj > spin-resolved 

t ^ 

densities and the total spin polarization, P = ^ t +"-i , 
[rf = J2y n<T (y)) as a function of magnetic field. The 
features in show a striking similarity with the fea- 
tures of the LDOS, and the behavior of the P follows 
that one of the DOS calculated in the Hartree approx- 
imation (c. f. Fig. [5] (a),(b)). This similarity is not 
coincidental. The regions with the enhanced LDOS cor- 
respond to compressible strips, and it is the compressible 
strips where the spin splitting of subbands takes place. 
Indeed, in the compressible region the subbands are only 
partially filled because Jfd < 1 there, and, therefore, the 
population of the spin-up and spin-down subbands can 
be easily changed. This population difference triggered 
by Zeeman splitting is enhanced by the Hubbard inter- 
action leading to different effective potentials for spin-up 
and spin-down electrons and eventually to the subband 
spin splitting. 

For a more detailed analysis let us follow an evolution 
of the band structure in Fig. |3ja)-(c). The right panel 
of Fig. [3] (a) shows a spin-resolved magnetoband struc- 
ture corresponding to the case when LL1 forms a com- 
pressible strip in the middle of the ribbon for the case of 
spinless electrons (c.f. the left panel of the figure). The 
Hubbard interaction pushes up the spin-down subbands 
above the window \E — Ep\ ~ 2fc^T such that it becomes 
depopulated, and the compressible strip in the middle is 
occupied by spin-up electrons only. As a result, the DOS 
at Ep of the spin-up electrons is enhanced, while that 
one of the spin-down electrons is zero, see Fig. 2] (d). 
All these lead to the difference in the electron densities 
and and the spin polarization in the ribbon as 
shown in Figs. H] (b) and (c). When the magnetic field 
increases such that the LL1 is pushed above the Ep, the 
compressible strip in the middle disappears, the DOS at 
Ep for both spin species becomes equal to zero and the 
spin polarization vanishes. This is illustrated in the band 
diagram shown in Fig. [3] (b) corresponding to the case 
when Ep is situated between LLO and LL1. With fur- 
ther increase of the magnetic field the LLO is pushed up 
to Ep. As this subband is pushed from below, in this case 
it is the higher energy spin-down state that gets pinned 
to Ep forming a compressible strips, whereas the spin-up 
subband in the middle of the ribbon remains below the 
\E — Ep\ ss 2fceT. As a result, for this case the DOS 
at Ep for the spin-down electrons is larger than that for 
the spin- up electrons, Fig. HJd). Note that despite of 
this, rv > nr (see Fig. |U(b)) because the spin-up sub- 
band is fully occupied, whereas the spin-down subband 
is occupied only partially. 

Figure 0|e) shows the effective g*-factor for the arm- 
chair ribbon defined according to g* = (V 7 — V^)/hbB, 
where the averaging is done over all carbon atoms. Be- 
cause the electron density is related to the potential, the 
features in the g*-factor resemble those of the polariza- 
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tion P, showing the behavior reflecting successive pop- 
ulation and depopulation of the spin-up and spin-down 
subbands. 

In high magnetic field, B > 50 T, (corresponding 
to population of LL1 and LLO) the effective g-factor 
varies between 2.1 < g* < 2.7 with the average value 
(g*) w 2.25, which represents a rather modest enhance- 
ment in comparison to the case of noninteracting elec- 
trons in pristine graphene with g = 2. Note that in the 
bulk graphene the effective <?*-factor was reported to be 
g* w 2.7J£ 

It is worth mentioning that the main features of the 
spin polarization and the subband evolution in magnetic 
field resemble those in the cleaved-edge overgrown quan- 
tum wires (CEOQW)^i. This is because in both cases 
the potential corresponds to the hard- wall confinement. 
The difference is that in CEOQW, as well as in con- 
ventional GaAs split-gate wires^, the polarization and 
thus the effective <?*-factor are enhanced by a factor of 
~ 10 in comparison to the Zeeman splitting, whereas in 
the armchair ribbons this enhancement it is just ~ 0.22. 
This can be explained by fact that the bare g-factor in 
armchair GNRs is much larger that that one in GaAs ( 
<?GaAs/<?graphene ~ 0.022), such that Zeeman interaction 
in graphene remains dominant in comparison to the ex- 
change one. 

One more important difference of the graphene ribbons 
from the conventional GaAs quantum wires is in the char- 
acter of spin polarized edge states in the vicinity of the 
boundaries. In the conventional quantum wires the edge 
state of opposite spins are spatially separated^ This is 
ultimately related to the formation of the compressible 
strips near the boundaries of the split-gate wire because 
of the soft confinement due to the gatesJ^ In GNR due to 
the hard- wall confinement the compressible strips do not 
form near the boundaries, and hence the spatial separa- 
tion of the edge states of opposite spins does not occur. 
It is worth mentioning that in this respect the GNRs are 
also similar to CEOQW.-^ 

Let us now turn to the case of zigzag nanoribbons. 
The main features of the spin polarization and the sub- 
band evolution are rather similar to the armchair ones. 
There is however one important difference related to the 
presence of the zero-energy mode residing at the zigzag 
edges. In contrast to other modes exhibiting successive 
population and depopulation of the Landau levels in the 
middle of the ribbon, this mode always stays pinned to 
Ep thus forming a compressible strip with the enhanced 
density of states at the edges for all magnetic fields . The 
Hubbard interaction leads to a complete spin polarization 
such that electrons in this mode are always in the spin-up 
state. This is seen in the spatially-resolved polarization 
shown in Fig. HJg). Because of this the total spin- up den- 
sity is significantly larger than the spin-down one for all 
magnetic fields, and the DOS for spin- up electrons never 
drops to zero. This is in contrast to the case of arm- 
chair GNRs where electron densities for opposite spin 
species can be equal when the Fermi energy lies between 



two consecutive Landau levels, and where both DOS for 
spin-ups and spin-downs can drop to zero (c.f. Figs, (d) 
with (j)). Because of the strong polarization near zigzag 
edges, the effective <7*-factor is strongly enhanced at low 
field (when the Zeeman spliting is small), and at higher 
fields it decreases to values g* w 3 comparable to those 
in the armchair GNRs, see Fig. |4](k). 

Figures|3If),(l) show evolution of the two-terminal con- 
ductance as the magnetic field increases for the arm- 
chair and zigzag ribbons. The conductance is appar- 
ently spin polarized with Q* ^ . This reflects the 
fact that at a given magnetic field the number of prop- 
agating states at Ep accommodating spin-up and spin- 
down electrons are different. As in the case of spinless 
electrons, the spin-resolved conductance also exhibits a 
bump-like structure which origin is the same as for the 
case of spinless electrons. Note that because of the pres- 
ence of a spin-polarized zero-energy mode, the spin-up 
conductance is larger than spin-down one for most fields. 
Finally we stress that Figs. EJf),(l) show the conduc- 
tances of ideal GNRs without defects. The defects scat- 
tering will modify GNRs conductance, especially for the 
low-velocity modes flowing in the middle of the ribbons. 
At the same time, the edge states corresponding to the 
classical skipping orbits are robust against the impurity 
scattering*^ 

Finally we note that we also performed similar com- 
putations of the spin polarization for GNRs using the 
density functional theory with the exchange functional 
proposed by Polini et a/4£ including the spin degree of 
freedom as prescribed in Ref. [4^ . Practically no spin 
polarization was observed that we attribute to the pos- 
itive sign of the exchange energy in Ref. l48l . More sys- 
tematic studies of the spin polarization in GNRs using 
different approach (such as the spin-DFT, Hartree-Fock, 
etc.) would be very interesting. 

Note that spin polarization and the enhancement of the 
g-factor in bulk graphene in the presence of impurities 
was recently studied by Volkov et al£H. 



IV. CONCLUSIONS 

We provide a systematic quantitative description of the 
spin polarization, the subband structure and the den- 
sity and potential profiles in the armchair and zigzag 
graphene nanoribbons in a perpendicular magnetic field. 
In our study we addressed realistically wide nanoribbons, 
and our conclusion concerning the density and potential 
distributions near the edge can also be applied for the 
case of a semi- infinite graphene sheet. Our calculations 
are based on the self-consistent Green's function tech- 
nique where electron interaction and spin effects are in- 
cluded by the Hartree and the Hubbard potentials. 

We first focus on the case of spinless electrons and find 
that the potential profile and the density distribution are 
different near the edges of the armchair and zigzag rib- 
bons. For the armchair termination, the potential at the 



8 



edge has a triangular shape whereas for the zigzag rib- 
bons it exhibits a well-type character. Both termina- 
tions show strong electron accumulation near the edges 
but this accumulation is stronger at the the zigzag edge. 
This difference is attributed to a topological property of 
the zigzag edge termination supporting the zero-energy 
mode. 

Because of the spin polarization in nanoribbons and 
conventional quantum wires is ultimately related to the 
formation of the compressible strips for the case of spin- 
less electrons, we study the LDOS, DOS, and the magne- 
tosubband structure for the armchair and zigzag ribbons 
focussing on the differences and similarities between them 
as well as on the formation of the compressible strips in 
the ribbons. For both types of nanoribbons we find that a 
compressible strip with the enhanced DOS forms in the 
middle of the ribbon in accordance with the successive 
population and depopulation of the Landau levels. For 
the case of zigzag edge termination we find a strong elec- 
tron accumulation and formation of a compressible strip 
near the edges over all the range of magnetic fields. This 
is caused by the presence of the zero-energy mode that 
always stays pinned to the Fermi energy because of the 
high DOS that this mode accommodates. 

Accounting for the Zeeman interaction and describing 
the spin effects via the Hubbard potential we discuss how 
the spin-resolved subband structure evolves when an ap- 
plied magnetic field varies. We find that the local spin po- 
larization of the electron density and the total spin polar- 
ization exhibit a behavior similar to that one of the LDOS 
and DOS for spinless electrons. This similarity is not co- 
incidental and reflects the fact that the regions with the 
enhanced DOS correspond to compressible strips where 



the spin splitting of subbands takes place. We find that 
for the armchair ribbons in high magnetic field the ef- 
fective g-factor varies between 2.1 < g* < 2.7 with the 
average value (<?*) « 2.25. For the zigzag nanoribbons 
we find that the zero-energy mode remains pinned to Ep 
and becomes fully spin-polarized for all magnetic fields, 
which, in turn, leads to a strong spin polarization of the 
electron density near the zigzag edge. Due to the contri- 
bution of the fully spin polarized zero-energy mode the 
effective <7*-factor in the armchair GNRs is strongly en- 
hanced at low field (when the Zeeman splitting is small) , 
and at higher fields it decreases to values g* « 3 compa- 
rable to those in the armchair GNRs. 

It is worth mentioning that the main features of the 
spin polarization and the subband evolution in magnetic 
field resemble those in the cleaved-cdgc overgrown quan- 
tum wires (CEOQW). This is because in both cases the 
potential corresponds to the hard- wall confinement. 

Finally, we stress the importance of accounting for 
the global electrostatics in the system at hand for the 
accurate description of the spin polarization in GNRs. 
(In the present study it is done by accounting for the 
long-range Coulomb interaction by means of the self- 
consistent Hartree potential). This is because the global 
electrostatic is responsible for the formation of the com- 
pressible strips, and it is the compressible strips where 
the spin splitting of subbands takes place. 
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